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ABSTRACT 

We have studied the atoll source 4U 1608-52 using a large data set obtained with the Rossi 
X-ray Timing Explorer. We find that the timing properties of 4U 1608-52 are almost exactly 
identical to those of the atoll sources 4U 0614+09 and 4U 1728-34 despite the fact that contrary 
to these sources 4U 1608-52 is a transient covering two orders of magnitude in luminosity. The 
frequencies of the variability components of these three sources follow a universal scheme when 
plotted versus the frequency of the upper kilohertz QPO, suggesting a very similar accretion flow 
configuration. If we plot the Z sources on this scheme only the lower kilohertz QPO and HBO 
follow identical relations. Using the mutual relations between the frequencies of the variability 
components we tested several models; the transition layer model, the sonic point beat frequency 
model, and the relativistic precession model. None of these models described the data satisfactory. 
Recently, it has been suggested that the atoll sources (among them 4U 1608-52) trace out similar 
three-branch patterns as the Z sources in the color-color diagram. We have studied the relation 
between the power spectral properties and the position of 4U 1608-52 in the color-color diagram 
and conclude that the timing behavior is not consistent with the idea that 4U 1608-52 traces 
out a three-branched Z shape in the color-color diagram along which the timing properties vary 
gradually, as Z sources do. 



Subject headings: accretion, accretion disks 
neutron stars: oscillations — X-rays: stars 

INTRODUCTION 



binaries: close — stars: individual (4U 1608-52) — stars: 



Most of the neutron star low-mass X-ray bi- 
naries can be divided into two classes, Z and 
atoll sources, based on the correlated behav- 
ior of their timing properties at low frequencies 
{v < 200 Hz) and their X-ray spectral properties 
(Hasinger & van der Klis 1989). Both classes show 
quasi-periodic oscillations (QPOs) with frequen- 
cies ranging from a few hundred Hz to more than 
1000 Hz (kilohertz QPOs). The low-frequency 
part of the power spectra of both classes is usu- 
ally dominated by a similar broad band-limited 
noise component, but it is unclear if these are 
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physically the same component (Wijnands & van 
der Klis 1999). In addition to the band-limited 
noise both classes show several quasi-periodic os- 
cillations below 200 Hz. In the Z sources these 
are named after the branch of the Z track in 
the color-color diagram (see below) where they 
mostly occur: horizontal (HBOs), normal (NBOs) 
and fiaring branch oscillations (FBOs). The HBO 
shows a sub- and a second harmonic at ~ 0.5 
and ~ 2 times the frequency of the main peak, 
and sometimes a peak at ^ 1.5 times the fre- 
quency of the HBO (Jonker et al. 2002). Obvi- 
ously, the most straightforward interpretation of 
this is that the sub HBO is the fundamental fre- 
quency, with second, third and fourth harmonics 
all observed (Jonker et al. 2002). Below 200 Hz 
the atoll sources show several Lorentzian compo- 
nents (see e.g. van Straaten et al. 2002). Note 
that these components are called Lorentzian and 
not QPO; this is because these features are some- 
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times too broad (FWHM > centroid frcqucncy/2) 
to be classified as a QPO. In the atoll sources all 
components become broader as their character- 
istic frequency decreases (Psaltis, Belloni & van 
der Klis 1999; van Straaten et al. 2002). It has 
been suggested that the HBO in the Z and the 
low frequency Lorentzian in the atoll sources are 
the same physical components (Ford & van der 
Klis 1998; Psaltis et al. 1999; Wijnands & van der 
Klis 1999). 

The energy spectrum of neutron star low- 
mass X-ray binaries can be usefully parametrized 
through the use of color color or color intensity 
diagrams, where a color is the ratio of counts in 
two different energy bands. As the energy spec- 
trum of a source changes, it moves through these 
diagrams. The timing properties of both the Z 
and the atoll sources depend on the position of the 
source in the color-color diagram. The Z sources 
move fast through the color-color diagram and 
draw up a Z track within hours to days (sec e.g. 
GX 340-1-0; Jonker et al. 2000) whereas the typical 
atoll sources move slowly through the color-color 
diagram and draw up a C-shaped track within 
weeks to months (sec e.g. 4U 1728-34; Di Salvo 
et al. 2001), although the "banana" part of the 
track, a curved branch to the bottom and right 
hand sides of the diagram, is traced out as fast as 
in the Z sources. The slow motion in the island 
part of the diagram, to the left hand and top sides, 
combined with observational windowing, tends to 
lead to the formation of isolated patches of data 
points, which is the origin of the term "island 
state" (Hasinger & van der Klis 1989). Whereas 
the continuum power spectra of the banana state 
are dominated by a power law component at low 
frequencies with perhaps a weak band-limited 
noise component (which becomes stronger as the 
source approaches the island state), the island 
state power spectra arc dominated by the band- 
limited noise. As the source moves away from 
the banana into the island state the count rate 
drops, the X ray spectrum becomes harder and 
the band-limited noise becomes stronger while its 
characteristic frequency decreases. In the most 
extreme island states (4U 1608-52; Yoshida et al. 
1993; Mendez et al. 1999; 4U 0614-^09; Mendez 
et al. 1997; van Straaten et al. 2000; 4U 1728-34; 
Ford & van der Klis 1998; Di Salvo et al. 2001; 4U 
1705-44; Langmeier, Hasinger & Triimper 1989; 



Bcrgcr & van der Klis 1998; Ford, van der Klis & 
Kaaret 1998; Aql X-1; Reig et al. 2000; KS 1731- 
260; Muno et al. 2000) the source is faint and hard 
and the band limited noise is strong; this is the 
state in which neutron stars are most similar to 
BHCs in the low hard state (van der Klis 1994a; 
Crary et al. 1996; Olive et al. 1998; Bcrgcr & van 
der Klis 1998). Because of the low count rates and 
slow motion through the color-color diagram the 
precise properties of the extreme island state have 
been hard to ascertain, although observations of 
4U 1608-52 with Tenma (Mitsuda et al. 1989) 
indicated the existence of an extended branch in 
this state that was traced out over an interval of 
weeks. 

Recently, Muno, Remillard & Chakrabarty 
(2002a) and GierUnski & Done (2002a) used large 
data sets from the Rossi X-ray Timing Explorer 
(RXTE) to study the color— color diagrams of sev- 
eral of the Z and atoll sources, which contain inter- 
esting additional information about, in particular, 
the nature of the extreme island state in the tran- 
sient atoll sources 4U 1608-52, 4U 1705-44 and 
Aql X-1. They suggested that the atoll sources 
trace out similar three-branch patterns as the Z 
sources, with the extreme island state cast in the 
role of Z source horizontal branch. However, these 
authors did not address the timing properties of 
the sources they studied. 

In this paper, we continue our work on the cor- 
related X-ray spectral and timing behavior of Z 
and atoll sources using RXTE (Wijnands et al. 
1996, 1997, 1998a,b; Ford et al. 1998; Mendez et 
al. 1997, 1999; Mendez & van der Klis 1999; Reig 
et al. 2000; .Jonker et al. 2000, 2002; van Straaten 
et al. 2000, 2001; Di Salvo et al. 2001; Homan et 
al. 2002) with an analysis of 4U 1608-52. In these 
previous analyses, we have found strong correla- 
tions between the behavior of the timing features 
and the position of the source in the color color 
diagram in both Z and atoll sources. 4U 1608- 
52 is a transient source that shows outbursts with 
a recurrence time varying between 80 days and 
several years (Loclmcr & Roussel-Dupre 1994). 
low-frequency QPOs were discovered in the "is- 
land" state of 4U 1608-52 with the Ginga satel- 
lite (Yoshida et al. 1993) and kilohertz QPOs were 
discovered with RXTE (van Paradijs et al. 1996; 
Berger et al. 1996). The source was included in 
the samples studied by Muno et al. (2002a) and 
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Gierlinski & Done (2002a) and was one of the 
sources that was reported to show Z-hke behav- 
ior in the color-color diagram. By connecting the 
timing with the energy spectral properties we can 
test whether 4U 1608-52 indeed behaves as a Z 
source. If this is true one would expect the power 
spectral properties to change smoothly along the 
Z track, as is the case in the Z sources. We find 
that this is not the case. In §5 we also perform a 
more general comparison of frequencies observed 
in the Z and the atoll sources. 

Parallel tracks in color-color and color-intensity 
diagrams were first observed in the Z sources 
(Hasinger et al. 1990; Kuulkers et al. 1994) and 
later also in the atoll sources 4U 1636-53 (Prins 
& van der Klis 1997; Di Salvo, Mendcz & van der 
Klis 2003), 4U 1735-44 (Wijnands et al. 1998c) 
and 4U 0614-h09 (van Straaten et al. 2000). Re- 
cently Muno et al. (2002a) reported further paral- 
lel tracks in the color-intensity diagrams of several 
atoll sources. The parallel-track phenomenon in 
the plot of lower kilohertz QPO frc;quency ver- 
sus intensity was first observed in 4U 0614-1-09 
(Ford et al. 1997) and 4U 1608-52 (Yu et al. 
1997) and has since been observed in several more 
atoll sources (see e.g. Mendez 2000). It has been 
suggested that the parallel tracks in the plot of 
intensity versus frequency of the lower kilohertz 
QPO and the parallel tracks in the color-intensity 
diagrams might be the same phenomena (van der 
Klis 2000, 2001; Muno et al. 2002a). Van der 
Klis (2001) has proposed a possible explanation 
for this parallel track phenomenon in terms of a 
filtered response of part of the X-ray emission 
to changes in the mass accretion rate through the 
disk. The frequency versus intensity parallel-track 
phenomenon is particularly clear in 4U 1608 52 
(Mendez et al. 1999). In §4 we investigate the 
relation between these frequency versus intensity 
parallel tracks and the parallel tracks in the color- 
intensity diagrams. 

2. OBSERVATIONS AND DATA ANAL- 
YSIS 

In this analysis we use all available public data 
from 1996 March 3 to 2000 May 24 from RXTE's 
proportional counter array (PCA; for more instru- 
ment information see Zhang et al. 1993). The data 
are divided into observations that consist of one to 



several satellite orbits. We exclude data for which 
the angle of the source above the Earth limb is 
less than 10 degrees or for which the pointing off- 
set is greater then 0.02 degrees. In our data set we 
found 7 type I X-ray bursts and we exclude those 
(~ 16 s before and ~ 240 s after the onset of each 
burst) from our analysis. 

We use the 16-s time-resolution Standard 2 
mode to calculate the colors. For each of the 
five PCA detectors (PCUs) we calculate a hard 
color, defined as the count rate in the energy band 
9.7—16.0 keV divided by the rate in the energy 
band 6.0—9.7 keV, and a soft color, defined as the 
count rate in the energy band 3.5—6.0 keV divided 
by the rate in the energy band 2.0—3.5 keV. Per 
detector we also calculate the intensity, the count 
rate in the energy band 2.0—16.0 keV. To obtain 
the count rates in these exact energy ranges we in- 
terpolate linearly between count rates in the PCU 
channels. We calculate the colors and intensity 
for each time interval of 16 s. We subtract the 
background contribution in each band using the 
standard background model for the PCA version 
2.1e. 

In order to correct for the changes in effective 
area between the different gain epochs and for the 
gain drifts within those epochs as well as the dif- 
ferences in effective area between the PCUs them- 
selves we used the method introduced by Kuulkers 
et al. (1994): for each PCU we calculate the colors 
of the Crab, which can be supposed to be constant 
in its colors, in the same manner as for 4U 1608- 
52. We then average the 16 s Crab colors and in- 
tensity per PCU for each day. In Figure 1 we can 
see the clear differences in the soft color trends of 
Crab between the five PCUs caused by the effects 
mentioned above. For each PCU we divide the 16 s 
color and intensity values obtained for 4U 1608 52 
by the corresponding Crab values that are closest 
in time but within the same gain epoch. We then 
average the colors and intensity over all PCUs. If 
we multiply the soft color by 2.36, the hard color 
by 0.56 and the intensity by 2400 c/s/PCU we 
approximately recover the observed, uncorrected 
values for 4U 1608 52 (quoted values are averages 
of the Crab colors during the observations over all 
five PCUs). To improve statistics we rebin the 
16 s color and intensity points to 256 s and ex- 
clude all data for which the resulting relative er- 
rors are larger than 5%. This led to a loss of about 
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40 ks of data at count rates below ~ 0.01 Crab. 
For most of these data neither timing nor colors 
were of sufficient quality to determine the source 
state. Starting May 12, 2000, the propane layer 
on PCUO, which functions as an anti-coincidence 
shield for charged particles, was lost. However, as 
all our data after May 12, 2000 (^10 ks) was ex- 
cluded because the relative color errors were larger 
than 5% this was not a issue in our analysis. 

For the Fourier timing analysis we use the 
122— /is time-resolution Event and Single Bit 
modes and the 0.95— /xs time-resolution Good 
Xenon modes. We only use data for which all 
energy channels are available. This led to a loss of 
about 9 ks of data in the lower banana state. In 
some observations there was data overflow in the 
timing modes due to excessive count rates from the 
source; we excluded those data. This led to a loss 
of about 27 ks of data at the highest intensity level 
(see also below). The power spectra were con- 
structed using data segments of 256 s and 1/8192 s 
time bins such that the lowest available frequency 
is 1/256 Hz and the Nyquist frequency 4096 Hz; 
the normalization of Leahy et al. (1983) was used. 
To get a first impression of the timing properties 
at the different dates and positions in the color- 
color diagram, the power spectra were combined 
per observation. The resulting power spectra were 
then converted to squared fractional rms. We sub- 
tracted a constant Poisson noise level estimated 
between 2000 and 4000 Hz where neither noise nor 
QPOs are known to be present. As fit function 
we use the multi-Lorentzian function; a sum of 
Lorentzian components (Grove et al. 1994; Olive 
et al. 1998; Belloni, Psaltis & van der Klis 2002) 
plus a power law to fit the VLFN (see van Straaten 
et al. 2002). We only include those Lorcntzians 
in the fit whose significance based on the error in 
the power integrated from to oo is above 3.0 u. 
Two to five Lorentzian components were generally 
needed for a good fit. For the Lorentzian that is 
used to fit the band-limited noise we fixed the cen- 
troid freqTiency to zero. We plot the power spectra 
in the power times frequency representation (yVi,; 
e.g. Belloni et al. 1997; Nowak 2000), where the 
power spectral density is multiplied by its Fourier 
frequency. For a multi-Lorentzian fit function 
this representation helps to visualize a character- 
istic frequency, i'max; namely, the frequency where 
each Lorentzian component contributes most of 



its variance per logarithmic frequency interval: in 
vVv the Lorentzian's maximum occurs at t'max 
(z^max = v^t'o + A^, where vq is the centroid and 
A the HWHM of the Lorentzian; Belloni et al. 
2002). We represent the Lorentzian relative width 
by Q defined as 

3. A TIMING TOUR THROUGH THE 
COLOR-COLOR DIAGRAM 

In this section we try to get an idea of the tim- 
ing properties of the source as it moves through 
the color-color diagram. We step through the 
data in chronological order and look at the tim- 
ing properties (per observation, see §2) and the 
position of the source in the color-color diagram. 
In this way we select continuous time intervals for 
which the power spectrum and the position in the 
color-color diagram remain similar. For each con- 
tinuous time interval we construct a representa- 
tive power spectrum by adding up all observations 
within that interval. To get a first idea of the 
timing properties we made initial fits to the rep- 
resentative power spectra of the continuous time 
intervals using the multi-Lorentzian fit function 
described in §2. In most cases the power spectral 
features present could be directly identified with 
known components seen in other atoll sources by 
comparing with the power spectral features of 4U 
0614+09 and 4U 1728-34 (using Figures 1, 2 and 3 
of van Straaten et al. 2002). A more thorough fit- 
ting and description of the power spectral features 
will be presented in §3.4 where we obtain optimal 
statistics by constructing representative intervals 
for all our data by adding up the continuous time 
intervals that have similar positions in the color- 
color diagram and show similar power spectra. 

The data on 4U 1608-52 can be usefully di- 
vided into 3 segments. The first segment ranges 
from 1996 March 3 to December 28 (the decay 
of the 1996 outburst, see Berger et al. 1996), the 
second segment from 1998 February 3 to Septem- 
ber 29 (the 1998 outburst, sec Mcndcz et al. 1998) 
and the third segment from 2000 March 6 to May 
10 (persistent data). In practice most data were 
available for the second segment (the 1998 out- 
burst) so we will present the results for the second 
segment first. These results can then serve as a 
template for the rest of the data. 
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3.1. Segment 2 (the 1998 outburst) 

In Figure 2 wc show as black points the 
lightcurve and color-color diagram for the sec- 
ond segment of the data (the 1998 outburst). The 
grey points in the lower frame represent the overall 
color-color diagram including all the data. Inten- 
sity and colors are normalized to the Crab (sec §2). 
The numbers 9 17 (1 8 arc reserved for the first 
segment of the lightcurve) indicate the continuous 
time intervals for which the power spectra and 
the position in the color-color diagram remain 
similar (see §3). The data for which no power 
spectra could be computed due to data overflow 
(see also §2) arc indicated by the larger crosses 
at the peak of the outburst and in interval 9 of 
the color-color diagram. In Table 1 we present 
the duration of each interval, the time until the 
next interval and the 2-16 keV intensity during 
each interval. We identify several power spectral 
features and classify the continuous time intervals 
into color intervals A- J; these will be discussed in 
§3.4 where we put all information together. For 
the 1998 outburst we find 7 different classes, i.e., 
the source returns twice to similar positions in 
the color-color diagram where it displays similar 
power spectral shapes. Each class is marked with 
a different symbol in the color-color diagram of 
Figure 2. 

For the 1998 outburst wc confirm the result of 
Mendez et al. (1999) that the color-color diagram 
shows the classical atoll shape (see Hasinger & van 
der Klis 1989). The characteristic frequencies of 
most of the power spectral components increase 
along the track starting at the open triangles and 
ending at the crosses in Figure 2 (for a more de- 
tailed discussion of the power spectral behavior 
with respect to the position in the color-color di- 
agram, see §3.5). 

3.2. Segment 1 (the decay of the 1996 out- 
burst) 

In Figure 3 wc show as black points the 
lightcurve and color-color diagram for the first 
segment of the data (the decay of the 1996 out- 
burst). The grey points represent the overall 
color-color diagram including all the data. The 
numbers 1-8 again indicate continuous time inter- 
vals for which the power spectra and the position 
in the color-color diagram remain similar. The 



symbols and entries in Table 1 are as described in 
§3.1. 

We find one class additional to those observed 
during intervals 9-17. This class is composed of 
continuous time intervals 2, 4 and 6, fills up the 
region between the open diamonds and the filled 
stars in the color-color diagram, and is marked 
with open stars in Figure 3. Note further that part 
of the open diamonds here are in a slightly differ- 
ent position in the color-color diagram, at a higher 
hard color and a lower soft color, than where they 
were during the second segment of the lightcurve. 
This is accompanied by a lower intensity. If we 
sort the classes by characteristic frequency as mea- 
sured in the power spectra, the characteristic fre- 
quencies change along an "e-shaped" track in the 
color-color diagram (see also §3.5). 

3.3. Segment 3 

In the third segment of the data the source 
count rates were low and in most observations 
one or more detectors were switched off. This, to- 
gether with the short durations of the observations 
(only ~ 500 — 2500 s per observation) led to bad 
statistics. In most cases it was impossible to iden- 
tify any power spectral features, therefore the clas- 
sification for this segment of the data was solely 
done on position in the color-color diagram. If any 
power spectral features were detected they were 
always consistent with this color-based classifica- 
tion. As before, in Figure 4 black points repre- 
sent the lightcurve and color color diagram for the 
third segment of the data and grey points the over- 
all color-color diagram. The numbers 18-30 now 
indicate continuous time intervals for which the 
position in the color-color diagram remained sim- 
ilar, and the symbols mark the classification, here 
based on the position in the color color diagram 
only; entries in Table 1 are as described in §3.1. 
Note that the source is in the class marked with 
the crosses for only the second time during con- 
tinuous time interval 21, but now at much lower 
intensity (0.03 Crab) than it was the first time in 
continuous time interval 9 (0.24-1.30 Crab) at the 
peak of the 1998 outburst. 

3.4. The combined power spectra 

To improve the statistics wc average the power 
spectra in each of the 8 classes. In the class marked 
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with the fiUed circles the lower kilohertz QPO is 
extremely narrow and varies in frequency over sev- 
eral hundreds of Hz. This leads to multiple peaks 
in the power spectrum. Therefore wc split up this 
class into three parts depending on lower kilohertz 
QPO frequency; the first with the lower kilohertz 
QPO ranging from 540 to 640 Hz, the second with 
the lower kilohertz QPO ranging from 640 to 710 
Hz and the third with the lower kilohertz QPO 
ranging from 710 to 900 Hz. A finer division would 
have compromised the statistics at low frequen- 
cies. The observations in this class where no lower 
(or upper) kilohertz QPO was detected were added 
based on position in the color-color diagram. In 
Figure 5 we show the resulting 10 intervals in the 
color color diagram and in Figure 6 wc show the 
corresponding hard color and soft color vs. in- 
tensity diagrams. Note that the data for which 
no power spectra could be computed due to data 
overflow (see §2 and 3.1) are not included in Fig- 
ures 5 and 6. We mark the different classes in 
order of increasing characteristic frequencies from 
A to J. We use letters here to avoid confusion with 
the numbered continuous time intervals of §3. 

We fit each interval with the multi-Lorentzian 
fit function described in §2. For the intervals 
where the kilohertz QPOs have sufficiently high 
frequencies not to interfere with the low -frequency 
features and vice versa, we fit the kilohertz QPOs 
between 500 and 2048 Hz and then fix the kilohertz 
QPO parameters when we fit the whole power 
spectra. This is for computational reasons only; 
the results are the same as those obtained with all 
parameters free. If the Q value of a Lorentzian 
becomes negative in the fit, we fix it to zero. No 
significantly negative Q's occured. The fits have a 
X^/dof below 1.4 for intervals A-G. For interval H 
the x^/dof is high (5.43); this is caused by the mo- 
tion of the lower kilohertz QPO from 710 to 900 
Hz (see above). For intervals I and J the x^/dof 
are 2.5 (dof = 100) and 1.8 (dof = 97). This is 
caused by deviations of the VLFN from the power 
law used to fit it (see appendix A). 

We now specify the terminology that we will use 
for the various power spectral components. Typ- 
ical power spectral components observed for the 
atoll sources are described for the two sources 4U 
0614-^09 and 4U 1728-34 in van Straaten et al. 
(2002). As in Bclloni et al. (2002) we call the 
upper kilohertz QPO L„ and its characteristic fre- 



quency Vu- The lower kilohertz QPO has been 
linked to a broad bump at ~ 10 to ~ 25 Hz found 
in the low luminosity bursters IE 1724-3045, GS 
1826 24 and SLX 1735 269 (Bclloni et al. 2002) 
and in the atoll source 4U 0614-1-09 at its low- 
est characteristic frequencies (van Straaten et al. 
2002). Although we emphasize that this identifi- 
cation is very tentative (see §5.2.1) we will, as was 
done in Belloni et al. (2002), call both the lower 
kilohertz QPO and this 10-25 Hz bump (which do 
not occur simultaneously), (characteristic fre- 
quency Ui). No standard terminology yet existed 
for the hcctohertz Lorentzian (Ford & van der Klis 
1998); here we will call it Ij^hz (characteristic fre- 
quency UhHz)- The behavior of the band-limited 
noise and QPOs in the 0.1 50 Hz range in the atoll 
sources 4U 1728-34 and 4U 0614-^09 is complex 
(see also Di Salvo et al. 2001; van Straaten et al. 
2002), wc describe it as a function of or position 
in the color-color diagram (see Figure 8). First, 
the characteristic frequency of the band-limited 
noise increases from ^ 0.3 to ^ 15 Hz. In this 
phase the band-limited noise is broad and usu- 
ally fitted with a zero-centered Lorentzian or a 
broken power law. Then at 15 Hz this noise 
component appears to "transform" into a narrow 
QPO (called very low-frequency Lorentzian in van 
Straaten et al. 2002) whose frequency smoothly 
continues to increase up to ~ 50 Hz, while what 
appears to be another band-limited noise compo- 
nent appears at lower characteristic frequencies. 
It is unclear exactly how these three components 
are related. Here, to remain true to the naming 
scheme of Bclloni ct al. (2002) wc call the band- 
limited noise component that becomes a QPO as 
well as the QPO it becomes Lb (characteristic 
frequency vj,) and the "new" broad band-limited 
noise appearing at lower frequency L(,2 (charac- 
teristic frequency Vb2), but please note the un- 
certainties in the interpretation underlying this 
nomenclature. Finally, in the banana branch of 
4U 1728-34 and 4U 0614-^09 (where usually no 
kilohertz QPOs are detected) a broad Lorentzian 
is also present at ~ 20 — 35 Hz for which it is un- 
clear whether it is L;,, L(,2 or a new component; 
here we shall list it as L;,2. The atoll sources 4U 
1728-34 and 4U 0614-^09 both additionally show 
a ~ 1 — 50 Hz Lorentzian at frequencies above 
Lft, called the low-frequency Lorentzian by van 
Straaten et al. (2002). At high characteristic fre- 
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qucncics (i^max > 25 Hz) this Lorcntzian appears 
as a narrow QPO. At low characteristic frequen- 
cies (fmax < 10 Hz) this Lorentzian is broad and 
can be identified with the component in the low lu- 
minosity bursters which Belloni et al. (2002) call 
Lfi ('hump'). For this reason, we also call this 
low frequency Lorentzian (characteristic fre- 
quency Uh). In some observations the low lumi- 
nosity bursters show a narrow low-frequency QPO 
(called Llf in Belloni et al. 2002) simultaneously 
with L/(. The fmax of Llf is slightly lower than 
that of Lh- For IE 1724-3045 the Lorentzian cen- 
troid frequencies of Llf and L/i coincided but this 
was not the case for GS 1826-24 (Belloni et al. 
2002). Llf was not observed in 4U 1728-34 or 
4U 0614+09 and we also do not detect it in 4U 
1608-52. 

A detailed description of the behaviour of the 
power spectral components of 4U 1608-52 is given 
in appendix A. In Figure 7 we show the averaged 
power spectra and best fit functions of the inter- 
vals A- J. The fit parameters for the low-frequency 
part of the power-spectra are listed in Table 2, 
the fit parameters for the high-frequency part are 
listed in Table 3. We obtain 95% confidence upper 
limits for the fractional rms amplitude of LhHz in 
intervals B, C and E and of Lh in intervals B, C 
and E using Ax^ = 2.71, fixing Q to 0.2 and al- 
lowing VhHz to run between 100 and 200 Hz. For 
setting upper limits to Lh we fix Q to 1.5, 3.0 
and 3.5 and let Vh run between 22-32 Hz, 34-44 
Hz and 42-50 Hz respectively (the values expected 
from the results for 4U 1728-34 and 4U 0614-1-09). 

In Figure 8 we plot the characteristic frequen- 
cies of 4U 1608-52 versus together with the 
resuhs of van Straatcn et al. (2002) for 4U 1728- 
34 and 4U 0614+09. The black points mark the 
results for 4U 1608 52, the grey points the results 
for 4U 1728 34 and 4U 0614+09. The different 
symbols indicate different power spectral compo- 
nents. The results for 4U 1608 52 mostly fall on 
the relations established for 4U 1728 34 and 4U 
0614+09. Also the Q values of the various compo- 
nents versus Vu are similar to those for 4U 1728- 
34 and 4U 0614+09. In Figure 9 we plot the rms 
fractional amplitude of all components versus Vu- 
The general trends in this plot for 4U 1608-52 are 
similar to those of 4U 1728-34 and 4U 0614+09 
(van Straaten et al. 2002), but there is an offset 
in rms between the relations for the three sources 



(Mcndez, van der Klis & Ford 2001). For the two 
kilohertz QPOs the rms fractional amplitudes in 
4U 1608-52 are similar to those in 4U 0614+09 
but larger than in 4U 1728 34 (Di Salvo et al. 
2001). For all other features; LhHz, Lh, Lf,, Lb2 
and at low frequencies L^ the rms fractional am- 
plitudes for 4U 0614+09 arc the largest followed by 
4U 1728-34 whereas the rms fractional amplitudes 
for 4U 1608-52 are the lowest. As we describe in 
appendix A in more detail, both the relations in 
Figures 8 and 9 as well as a direct comparison with 
the power spectra of 4U 1728-34 and 4U 0614+09 
(Figures 1 and 2 in van Straaten et al. 2002) al- 
low us to identify the power spectral components 
of 4U 1608-52 in the different intervals within the 
identification scheme described above. 

3.5. How does 4U 1608-52 move through 
the color— color diagram? 

Now that we have described and identified all 
power spectral components in 4U 1608-52, we can 
link the timing properties (Tables 2 and 3) to 
the position in the color- color diagram. We can 
also link the timing properties and X-ray spectral 
properties to the classical island and banana states 
described by Hasinger & van der Klis (1989). In- 
tervals A-D are occurrences of the island state 
(strong broadband noise with low i^max, no VLFN) 
in which intervals A-C can be classified as the ex- 
treme island state. The extreme island state shows 
Lb, Lft, L^ and L„ all at low characteristic frequen- 
cies. The power spectral components are all broad 
and strong. Interval E forms a transition between 
the island state and the lower banana, it still has 
relatively strong broad band noise but also shows 
a pair of narrow kilohertz QPOs that are typical 
for what Hasinger & van der Klis (1989) called the 
lower left banana state. So, intervals F-H are all 
lower left banana state. Here the VLFN appears, 
the kilohertz QPOs are double, L(, transforms from 
a band-limited noise component into a QPO and 
L&2 appears (see also §3.4). Intervals I and most 
of interval J are also still in what Hasinger & van 
der Klis (1989) called the lower banana state based 
on their position in the color-color diagram; only 
the upper right part in the color-color diagram of 
interval J is in the upper banana state based on 
color-color position. The VLFN is strong and the 
broad band noise is weak in these intervals. Inter- 
vals I and J show L„ at the highest frequencies. 
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Lt,2, and a strong VLFN component. The charac- 
teristic frequencies increase in order A-J and form 
an "e-shaped" track in the color-color diagram. 

Recently, Muno et al. (2002a) and Gierlinski 
& Done (2002a) studied the color-color diagrams 
of several of the Z and atoll sources, including 4U 
1608 52, and suggcstc;d that the atoll sources trace 
out similar three-branch patterns as the Z sources. 
We observe the same shape of the color-color di- 
agram for 4U 1608-52 as Muno et al. (2002a) and 
Gierhnski & Done (2002a) did (Fig. 5). Inter- 
val C in Figure 5 represents a deviation from the 
classical atoll shape. According to the interpre- 
tation of Muno et al. (2002a) and Gierlinski & 
Done (2002a) the source would have to move in 
the Z-track order D-AB C in Figure 5. Interval 
C would then correspond to the horizontal branch 
of the Z sources. We observe a transition from D 
to C and back to D with gaps of three days dur- 
ing the decay of the 1996 outburst (Fig. 3). We 
also observe a transition from E to C and back to 
E with gaps of only one day during the decay of 
the 1998 outburst (Fig. 2). The source is first in 
E for an interval of ^ 0.1 day during which the 
hard color increases by ^ 0.04 (so Whard color ~ 0.4 
day~^), then there is a gap of 0.9 day after which 
the source appears in C with an increase in hard 
color of ^ 0.3 (if the source moved directly from E 
to C 

■fhard color ~ 0.3 day ^). Then in C there is 
an interval of ~ 0.16 day where the hard color in- 
creases by ~ 0.04 (fhard color « 0.3 day~^). So this 
is consistent with the source moving directly from 
E to C at approximately constant speed (proba- 
bly through D) and not through A or B as re- 
quired if 4U 1608-52 behaved as a Z source. In 
C there is a gap of ~ 1.8 day after which the 
source appears again in C but at a slightly 0.03) 
lower hard color. For an ~ 0.1 day interval the 
source stays in C while the hard color decreases 
by ^ 0.03 (uhard color ~ —0.3 day^"^). Then there 
is a gap of ~ 1.0 day and after that the source 
has returned to E with a hard color of ~ 0.3 less 
(I'hard color ~ —0.3 clav^''"). This is consistent with 
the source moving directly back from C to E, again 
at constant speed, and not through A or B. Tran- 
sitions from and to states A and B had gaps of 7 
and more days, so it is impossible to draw conclu- 
sions from these. The characteristic frequencies of 
the timing features decrease in the order D-C-B- 
A (appendix A). This is also not consistent with 



the idea that 4U 1608-52 behaves as a Z source, 
as in Z sources the characteristic frequencies of 
the timing features increase along the Z starting 
at the horizontal branch (i.e., this would predict 
frequencies decreasing in the order DBAC). One 
might say that 4U 1608-52 draws up an approx- 
imate "e" shape in the color color diagram along 
which the characteristic frequencies of the timing 
features change smoothly. 

4. A DETAILED STUDY OF THE PAR- 
ALLEL TRACK PHENOMENA 

Parallel tracks in the intensity versus lower kilo- 
hertz QPO frequency in 4U 1608-52 were discov- 
ered by Yu et al. (1997), and extensively studied 
by Mendez et al. (1999) and Mendez et al. (2001). 
The hard color and soft color vs. intensity dia- 
grams of 4U 1608-52 (see Fig. 6) and those of 
several other atoll sources also show narrow par- 
allel tracks (Prins & van der Klis 1997; Di Salvo 
et al. 2003; Muno et al. 2002a). Based on the ap- 
pearance of the diagrams, it has been suggested 
that there is a relation between these two types of 
parallel tracks (van der Klis 2000, 2001; Muno et 
al. 2002a). 

Based on our new analysis it is now possible to 

directly link the parallel tracks in the color inten- 
sity diagrams to those in the QPO frequency- 
intensity diagrams using the frequency of the 
lower kilohertz QPO as obtained by Mendez et 
al. (2001). Similarly to Mendez et al. (1999) and 
Mendez et al. (2001) we only include data for 
which both kilohertz QPOs are detected simulta- 
neously and therefore the lower kilohertz QPO is 
identified unambiguously. We rebin our 16 s colors 
in such a way that they match the 64-448 s data 
intervals of Mendez et al. (2001). In Figure 10 
we plot the frequency of the lower kilohertz QPO 
and the hard and soft color versus the 2.0-16.0 
keV intensity. The alternating black/grey sym- 
bols represent the parallel tracks in the intensity 
versus the frequency diagram; the tracks contain 
data that are continuous in time, with only the 
^ 2500 s gaps due to Earth occultations. Note 
that the parallel track phenomenon in the lower 
kilohertz QPO frequency vs. intensity diagram is 
only observed in a small region of the color in- 
tensity diagrams (see Figiire 6) where the lower 
kilohertz QPO is strong and narrow enough to be 
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accurately traced on 64-448 s timescales. 

With this analysis, we can investigate how 
the parallel tracks in the frequency vs. inten- 
sity diagram relate to those in the color inten- 
sity diagrams. Interestingly, it turns out that 
the frequency-intensity tracks can not, as previ- 
ously thought, be identified with the narrow ver- 
tical tracks visible in the color intensity diagrams 
of Figures 6 and 10. Each of these latter tracks 
corresponds to a single satellite orbit. The paral- 
lel tracks in the frequency vs. intensity diagram 
are only identifiable as (rather fuzzy) tracks in the 
color intensity diagrams after several satellite or- 
bits have elapsed and thus turn out to be com- 
posed of several of those narrow vertical tracks. 
The parallel tracks in the frequency versus in- 
tensity diagram could also identified in the hard 
color-intensity diagram of 4U 1636-53, where two 
banana tracks shifted by 20 % in intensity were 
observed in the hard color-intensity diagram (Di 
Salvo et al. 2003), and exactly the same 20 % shift 
was observed between the corresponding two par- 
allel lines in the frequency versus intensity dia- 
gram. In 4U 1608-52 we do not observe complete 
banana tracks shifted in intensity as in 4U 1636- 
53. The drift in X-ray flux in 4U 1608-52 already 
occurs on timescales of several days (Mendez et 
al. 1999), too short for a complete track to form, 
where in 4U 1636-53 the source can stay on one 
track for several months (Di Salvo et al. 2003). 
These fast changes in intensity in 4U 1608-52 are 
due to its transient character, where 4U 1636-53 
is a persistent source for which changes in inten- 
sity occur much more gradually. Note, that within 
each of the two parallel lines in the frequency ver- 
sus intensity diagram of 4U 1636-53 there is no 
correlation observed between frequency and inten- 
sity. 

The narrow vertical tracks in the color intensity 
diagrams are caused by variations in color that for 
the most part can be explained by the scatter due 
to counting statistics. The parallel tracks in the 
frequency versus intensity diagram, although com- 
posed of data covering several consecutive satellite 
orbits, already show up within individual orbits 
due to a short term correlation between the lower 
kilohertz QPO frequency and intensity in combi- 
nation with small errors in these quantities but 
persist over several consecutive orbits. The cor- 
relation between colors and intensity (or kilohertz 



QPO frequency) in individual orbits is veiled by 
the limited counting statistics, so in the color in- 
tensity diagrams the parallel QPO tracks can only 
be identified on the longer timescales of several 
satellite orbits. 

To further illustrate this point in Figure 11 we 
plot soft color versus the frequency of the lower 
kilohertz QPO; lines connect the points of indi- 
vidual satellite orbits, and for clarity four repre- 
sentative individual orbits are highlighted in black. 
On orbit time scales, the lower kilohertz QPO fre- 
quency and soft color seem uncorrelated due to 
the counting statistics errors on the color, whereas 
on a longer time scale a clear correlation emerges. 
We note that for some parallel tracks in Figure 
6 there is a correlation between color and inten- 
sity on short timescales, especially in interval J. In 
those cases the counting statistics are better due to 
the higher count rates and thus the parallel track 
phenomenon in the color diagrams is more simi- 
lar to that in the kilohertz QPO frequency versus 
intensity diagram (Muno et al. 2002a). However, 
in our interval J there is no timing feature present 
that is trackable on short timescales. 

In conclusion we can now identify the parallel 
QPO tracks in the color intensity diagrams. Bad 
counting statistics that shows up in the form of the 
narrow nearly vertical parallel tracks in the color 
intensity diagrams cause the color correlation 
to be veiled on short timescales (less than hours) 
so in the color intensity diagrams the parallel QPO 
tracks can only be found for the longer timescales 
of several satellite orbits. 

5. DISCUSSION 

We have studied the color diagrams and the 
power spectral behavior of the atoll source 4U 
1608-52. We found that the timing behavior of 
4U 1608 52 is almost identical to that of the other 
atoU sources 4U 1728-34 and 4U 0614-^09. If we 
plot the characteristic frequencies of the timing 
features versus the characteristic freqiiency of the 
upper kilohertz QPO, together with the results of 
van Straaten et al. (2002) for 4U 1728-34 and 4U 
0614+09, the three sources follow the same re- 
lations (see Fig. 8). Also the behavior of the 
Q value is the same for the three sources. The 
general trends in rms fractional amplitude for 4U 
1608-52 are also similar to those of 4U 1728-34 
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and 4U 0614+09, but there is an offset between 
the relations for the three sources (see Fig. 9). 
We connected the timing behavior with the posi- 
tion of the source in the color color diagram and 
found that the timing behavior is not consistent 
with the idea that 4U 1608-52 traces out a three- 
branched Z shape in the color-color diagram along 
which the timing properties vary gradually as is 
the case in in Z sources. Instead, the power spec- 
tral properties change along an "e shaped" track. 
Finally, our measurements for the colors together 
with the precise measurements of lower kilohertz 
QPO frequency of Mcndcz ct al. (2001) gave us 
an opportunity to link the parallel tracks in the 
intensity versus color diagrams with the paral- 
lel tracks in the intensity versus lower kilohertz 
QPO frequency diagrams. We found that the par- 
allel tracks in the frequency versus intensity di- 
agram can be found back in the color-intensity 
diagrams as fuzzy structures; they should not be 
confused with the narrow vertical parallel tracks 
visible in the intensity versus color diagrams which 
are mostly due to the errors in the colors. 

5.1. Comparison with other sources; col- 
ors 

The power spectral properties in 4U 1608-52 
change along an "e-shaped" track. The question 
now arises whether the "e-shaped" track we find 
in the color-color diagram of 4U 1608-52 is univer- 
sal for atoll sources. We can compare our results 
with those of Olive, Barret & Gicrlinski (2003) 
who studied the X-ray color and timing properties 
of a well sampled state transition of 4U 1705-44 
where the source moves from the lower banana to 
the extreme island state and back. This transition 
was also included in Muno et al. (2002a) who stud- 
ied a larger dataset of 4U 1705-44. In the color- 
color diagram the source moves from the lower ba- 
nana, to the left of the extreme island state as the 
count rate decreases. Then both the soft color and 
the count rate increase, whereas the hard color 
remains approximately constant, and the source 
traces out a horizontal track in the color-color di- 
agram. Then, while the count rate continues to in- 
crease, the source moves back from the right of the 
extreme island state to the lower banana (Muno et 
al. 2002a; Barret & Olive 2002; Olive et al. 2003). 
In the whole horizontal track the power spectra 
remained the same (Olive et al. 2003). 



In this extreme island state 4U 1705-44 did not 
reach characteristic frequencies as low as those in 
interval A of 4U 1608-52 (it is very similar to in- 
terval B). If we look at the overall light and color 
curves of 4U 1705-44 presented in Muno et al. 
(2002a) we see that apart from the extreme is- 
land state just described, which took place in the 
forty days around MJD 51230, there is another ex- 
treme island state, with much sparser sampling, in 
which a higher hard color is reached. If we take a 
quick look at the timing of this MJD 51380 data we 
find that when the source reaches this higher hard 
color, the source shows a power spectrum very sim- 
ilar to that of interval A of 4U 1608-52. As for 
the MJD 51230 extreme island state, during this 
MJD 51380 extreme island state the count rate 
and the soft color increase simultaneously and in 
the color-color diagram another horizontal track 
is drawn up above the MJD 51230 one. 

Returning now to 4U 1608-52 we note that in 
the state transition of 4U 1608-52 during the 1998 
outburst (see Figure 2) we observe the same phe- 
nomenon of a horizontal track being traced out in 
the extreme island state. The source is in interval 
C (continuous time interval 12 in Fig. 2) for two 
observations, the first of which has a count rate 
about a factor two higher than the second. This 
higher count rate is accompanied by a 5 % higher 
soft color in the first observation. The power spec- 
trum is the same in both observations. So, like the 
case of 4U 1705-44, the change in count rate is 
accompanied by a correlated change in soft color 
which traces out a small almost (the hard color 
changes slightly, see above) horizontal track in the 
color-color diagram (continuous time interval 12 
in Fig. 2). 

A possible explanation for this behavior of atoll 
sources in the extreme island state is the so- called 
"secular motion" in the color color diagram. This 
phenomenon was first observed in Z sources, in 
which the Z shaped track in the color-color and 
color-intensity diagrams is traced out within sev- 
eral hours up to a day. On longer timescales the 
whole Z track shifts both in soft color and count 
rate (e.g. Hasinger et al. 1990; Kuulkers et al. 1994; 
Jonkcr et al. 2000, 2002; Homan et al. 2002; Muno 
et al. 2002a). The timing properties remain mostly 
unaffected by these shifts and are primarily deter- 
mined by the position along the Z track (e.g. Ku- 
ulkers et al. 1994; Jonker et al. 2000, 2002; Homan 
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et al. 2002). The same phenomenon has been ob- 
served in the banana state of the atoll source 4U 
1636-53 (Prins fc van der Klis 1997; Di Salvo et 
al. 2003). The horizontal tracks in the extreme 
island state of the atoll sources may be entirely 
caused by a secular motion, similar to that in the 
Z sources. Because the sources remain in one par- 
ticular island state (similar timing properties and 
hard color) for a long time (weeks to months), the 
slow process of secular motion has time to draw up 
a horizontal branch. In this horizontal branch the 
timing remains similar, as is the case for Z source 
secular motion. As noted by Olive et al. (2003), 
this aspect of the behavior can be explained by 
the scenario that van der Klis (2001) proposed to 
explain the parallel track phenomenon in the in- 
tensity versus lower kilohertz QPO frequency dia- 
gram. 

It seems that the Z shape in the color-color di- 
agram of some of the atoll sources is caused by 
transitions between the banana and the extreme 
island states that occm- at different soft color. Fig- 
ure 12 provides a schematic of what in our inter- 
pretation occurs in these sources. Several different 
horizontally extended extreme island branches ap- 
pear above each other in the color-color diagram 
at different hard color values. These are traced out 
at different epochs. To first order only the hard 
color determines the timing properties. Within 
each extreme island state the soft color changes 
in correlation with intensity. During state transi- 
tions the extreme island state is entered or left at a 
soft color value depending on intensity. A possible 
explanation for the fact that the timing remains 
similar while the intensity (and thus soft color) 
changes and that the state transitions can occur 
at different intensities (soft colors) is the scenario 
of van der Klis (2001) where the truncation radius 
of the disk, which determines the timing properties 
and therefore the state, is not set by the accretion 
rate through the disk but by the ratio of this ac- 
cretion rate over its long-term average (c.f. Olive 
et al. 2003). So, if a source shows a state tran- 
sition to an extreme island state at a minimum 
intensity which increases after that (as happened 
in the MJD 51230 state transition of 4U 1705 44 
Muno et al. 2002a; Barret & Olive 2002; Olive et 
al. 2003), the source will enter the extreme island 
state at a lower soft color and move to a higher soft 
color as the intensity increases. If a state transi- 



tion to an extreme island state takes place during 
a decay after an outburst of a transient source, we 
would expect that after the extreme island state is 
entered at a particular soft color, the source then 
moves to a lower soft color as it fades (this and the 
reverse state transition from an extreme island to 
the banana state during the outburst rise were ob- 
served in Aql X-1, see Reig, van Straaten & van 
der Klis 2003). In our interpretation, the shape we 
observe in the color- color diagram of 4U 1608 52 
is caused by observations of the source in extreme 
island states at different intensities and therefore 
at different soft colors. Extreme island states A 
and B are observed at higher intensities than ex- 
treme island state C, causing C to have a lower 
soft color. Contrary to the case in 4U 1705 44 
and Aql X-1 we do not observe a large change in 
soft color (or intensity) within each extreme island 
state and it is unclear whether this is because the 
source was not observed sufficiently long or dense 
enough in these extreme island states, or whether 
this is a property of 4U 1608-52. 

5.2. Comparison with other sources; tim- 
ing 

5.2.1. Atoll sources, low luminosity bursters, and 

black hole candidates 

The three atoU sources, 4U 1608-52, 4U 1608- 
52 and 4U 0614-1-09, we have studied up to now 
(van Straaten et al. 2002, this paper) show a sim- 
ilar and very distinct timing behaviour (Figures 
8 and 9). 4U 1608-52 differs from 4U 0614-h09 
and 4U 1728-34 in that 4U 1608 52 is a transient 
source whose intensity changes by about 2 orders 
of magnitude where this is only about 1 order of 
magnitude for 4U 0614+09 and 4U 1728-34. Al- 
though 4U 1608-52 has a much larger range in 
intensity (~ 0.006 - 1.0 Crab in the 2.0-16.0 keV 
range) than 4U 0614+09 (- 0.017-0.1 Crab), 4U 
1608-52 and 4U 0614+09 show a similar range in 
characteristic frequencies for the different power 
spectral components. For 4U 1608-52 ranges 
from 216 to 1061 Hz (excluding intervals I and J 
for which the identification of z/„ is uncertain see 
appendix A). This range in v-u is reached with a 
minimum intensity range of ~ 0.05 — 0.1 Crab (see 
Fig. 6). For 4U 0614+09 a similar i/„ range of 233 
to 1067 Hz is reached with a minimum intensity 
range of ~ 0.03—0.05 Crab. So, both sources reach 
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similar frequency ranges within a similar change in 
intensity. 4U 1728-34 has not yet shown a similar 
range in characteristic frequencies as 4U 0614+09 
and 4U 1608 52 (imtil now i/„ > 400 Hz; Di Salvo 
et al. 2001; van Straaten et al. 2002). Also no 2.0- 
16.0 keV intensities were published for the inter- 
vals used in Di Salvo et al. (2001) and van Straaten 
et al. (2002). 

There are also some differences in the timing 
behaviour; in particular in Lf, when it appears as 
a QPO. As a function of Vu, i^b is slightly higher in 
4U 1608-52 than in 4U 1728-34 and 4U 0614+09 
but it covers the same; range. The Q values for L(, 
here are similar for the three sources. The relation 
between Uf, and Vu in Figure 8 has a turnover at 
30 Hz in 4U 0614+09 (van Straaten et al. 2000); 
there might be a turnover at 45 Hz in 4U 1728-34, 
but in 4U 1608-52 no turnover has been observed. 
The plot of rms fractional amplitude of all compo- 
nents versus Vu (Fig. 9) shows an offset between 
the relations for the three atoll sources 4U 1608- 
52, 4U 0614+09, and 4U 1728-34. This offset is 
different for the low-frequency features, L^hz, ^h, 
Lb, Lb2, and at low frequencies Lg, than for the 
two kilohertz QPOs (see also §3.4). It seems that 
there are two groups, one containing the high and 
the other the low-frequency features. The differ- 
ence between the two groups might be due to ge- 
ometrical effects. The different offsets in Figure 9 
hint towards the Lorentzian at ~ 200 Hz, found in 
the lowest frequency states of 4U 1608-52 (interval 
A) and 4U 0614+09, being L„ (see appendix A), 
whereas it could have been identified as either L„ 
or LhHz based on its frequency (see appendix A). 
The identification of the ~ 200 Hz Lorentzian is 
important because the power spectra of 4U 1608- 
52 and 4U 0614+09 at the lowest inferred mass 
accretion rate closely resemble those of the low 
luminosity bursters IE 1724-3045, GS 1826-24, 
and SLX 1735 269 (BcUoni et al. 2002) and the 
millisecond X-ray pulsar SAX J1808.4-3658 (Wi- 
jnands & van der Klis 1998). These sources all 
seem to be atoll sources that arc only observed 
at low mass accretion rate. If the ~ 200 Hz 
Lorentzian is the upper kilohertz QPO then the 
high-frequency peaks in those sources are prob- 
ably also upper kilohertz QPOs. Note that al- 
though in 4U 1608-52 Lh already becomes unde- 
tectable at ^ 20 Hz whereas this happens at ~ 45 
Hz in 4U 1728-34 and 4U 0614+09, this is not 



a significant difference as the upper limits on the 
fractional rms of in 4U 1608-52 are consistent 
with this component still being present at those 
frequencies. 

The broad L^ component found in intervals 
A and B was previously found in 4U 1608-52 
with Ginga by Yoshida et al. (1993) and has also 
been identified in the atoll source 4U 0614+09 
(van Straaten et al. 2002). This component has 
also been found in the low-luminosity bursters IE 
1724-3045, GS 1826-24, and SLX 1735-269 (Bel- 
loni et al. 2002). The Lorentzian has been ten- 
tatively identified with the lower kilohertz QPO 
based on extrapolations of frequency-frequency 
relations (see Psaltis et al. 1999; Belloni et al. 
2002; van Straaten et al. 2002). The rms frac- 
tional amplitude of the broad Li is much higher 
than that of the lower kilohertz QPO (see Fig. 9), 
howevcir. a similar increase in rms is visible in that 
Figure for L„. But note that L„ is detected in all 
intervals where L^ is present as the lower kilohertz 
QPO in intervals H-E then dissapears in intervals 
D and C and only shows up again as the broad 
Lorentzian in intervals B and A. 

Interval B of 4U 1608-52 shows a narrow QPO 
with a characteristic frequency (2.458 Hz) between 
those of Lfc and Lh- This QPO is probably the 
same as the one discovered in the island state of 
4U 1608-52 by Yoshida et al. (1993): their power 
spectrum of the 1989 August 25-26 interval is very 
similar to our interval B, it has Li, with a char- 
acteristic frequency of ~ 1 Hz, L/j at ~ 6 Hz 
and a narrow QPO at ^ 2 Hz (sec Figure 4 of 
Yoshida et al. 1993). The low-luminosity bursters 
IE 1724 3045, GS 1826 24 (BcUoni et al. 2002) 
as well as the BHCs GX 339 4 (Nowak, Wilms 
& Dove 2002) and Cyg X 1 (Pottschmidt et al. 
2002) also show narrow QPOs with a character- 
istic frequency between those of Lh and that of 
Lh- For IE 1724 3045 the centroid frequency of 
the narrow QPO coincides with the centroid fre- 
quency of Lh, but for 4U 1608 52 this is not the 
case; here the centroid frequency of the narrow 
QPO (2.45 ± 0.04 Hz) is close to half the cen- 
troid frequency of Lh (4.58 ± 0.28 Hz). No similar 
narrow Lorentzians were fitted by van Straaten 
et al. (2002) for 4U 1728-34 or 4U 0614+09, but 
4U 1728-34 does show narrow residuals between 
the characteristic frequency of L;, and that of Lh- 
In Figure 13 we plot Vh versus the Vmax of these 
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narrow QPOs. For Cyg X 1 many narrow QPOs 
were fitted by Pottschmidt et al. (2002); we only 
plot those QPOs that have a characteristic fre- 
quency between i/b and i/h- The points of the low- 
luminosity bursters and most of those of GX 339-4 
line up in Figure 13. These QPOs were all labeled 
as Llf by BcUoni ct al. (2002). The two points of 
GX 339-4 that deviate from the line are from two 
observations where GX 339-4 showed two simul- 
taneous QPOs; the highest frequency QPO falls 
on the relation in Figure 13 while the lowest fre- 
quency QPO falls below it. If we fit the points of 
the low-luminosity bursters and GX 339-4 with a 
power law our result as well as those of Yoshida 
et al. (1993) for 4U 1608-52 fall below this rela- 
tion. So we can not identify this narrow QPO in 
4U 1608-52 with Llf- In the two observations 
where GX 339-4 showed two narrow simultaneous 
QPOs, the characteristic frequencies of the lowest 
frequency QPO as well as the results for Cyg X-1 
also fall well below the power law. It might be that 
these QPOs and those of 4U 1608 52 are related. 
These points fall close to a line that indicates half 
the low-luminosity burster relation which is also 
plotted in Figure 13. Note that if we use ccntroid 
frequency instead of fmax the relations in Figure 
13 worsen. 

5.2.2. Z sources 

In addition to comparing the color-color di- 
agram tracks and the associated power spectra, 
we can also compare the timing properties of Z 
sources with those of atoll sources by making sim- 
ilar plots of fmax of the different power spec- 
tral components versus v-^ as we did for the atoll 
sources 4U 1608-52, 4U 0614-^09 and, 4U 1728- 
34 in Figure 8. In Figure 14 we do this for the 
Z sources GX 5-1 (Jonker et al. 2002), GX 340+0 
(Jonker et al. 2000), GX 17-^2 (Homan et al. 2002) 
and Cyg X-2 (Kuznetsov 2002). We include the 
lower kilohertz QPO, the Low-Frequency Noise 
(LFN), the Horizontal Branch Oscillations (HBO), 
and the harmonic and sub-harmonic of the HBO, 
and plot these versus Vu- The grey symbols are the 
results of the atoll sources also displayed in Fig- 
ure 8, the black symbols represent the Z sources. 
Note that the broad-band noise in these Z sources 
was not fitted with a zero-centered Lorentzian but 
with a cutoff power law, P{u) oc u~"e~'^^'^'"^\ 
for GX 5-1, GX 340-1-0, and GX 17-1-2 and with 



a smooth broken power law, P(i^) cx i^~"[l -|- 
(i//i/b)^]~^, for Cyg X-2. This leads to expres- 
sions for fmax, defined as the frequency of maxi- 
mum power density in j^Pj/, of (1 — a)i'cut for the 
cutoff power law, and [(1 — a)/ {13 -\- a — l)]^/^;/;, 
for the smooth broken power law. 

From Figure 14 we confirm the identification 
of the HBO in the Z sources with L^^ in the atoll 
sources made by Psaltis et al. (1999). However, 
the suggestions that cither the LFN (van dcr Klis 
1994b; Jonker et al. 2000) or the sub-HBO (Wi- 
jnands fc van der Klis 1999; Jonker et al. 2000) 
in the Z sources might be similar to the classical 
broad-band noise in the island state of the atoll 
sources seems not to be supported by the relation 
of the characteristic frequencies of these compo- 
nents with Vu. At Vu ~ 750 Hz the band-limited 
noise component in the L(,-L„ relation is replaced 
by a QPO (see §3.4) and the band-limited noise 
component 'Lb2 appears and follows a new relation 
with Vu. The LFN points fall below the Lf, — L„ 
relation but seem to line up with Lf,2. The sub- 
HBO points fall above the L;, — L„ relation until 
Vu « 750 Hz where the QPO takes over from the 
band-limited noise component (see above). 

Based on this comparison of frequencies, the 
LFN might still be associated with Lf,2 and the 
sub-HBO might be related with Lf, when it is a 
QPO, although for 4U 0614-1-09 and 4U 1728-34 
no harmonic relation was present between L/j and 
Lb when it is a QPO (van Straaten et al. 2002). 
Note also that the LFN in GX 17+2 follows a dif- 
ferent relation compared to that in the other Z 
sources. Other differences between GX 17+2 and 
the other Z sources in Figure 14 arc that the LFN 
in GX 17+2 is peaked, the harmonic of the HBO 
is relatively strong and it shows a flaring branch 
oscillation (FBO) whereas the other sources show 
a flat LFN, a relatively weak harmonic of the HBO 
and no FBO (Jonker et al. 2002). We note that 
none of these differences with the other Z sources 
serve to make GX 17+2 more similar to the other 
atoll sources. 

5.2.3. Summary 

The results presented in sections 3.4, appendix 
A, and this section show that the low-frequency 
part of the power spectra of both the atoll and the 
Z sources behaves in a very complex manner. For 
the atoll sources many different Lorentzian com- 
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ponents appear, disappear or change from broad 
Lorentzians into narrow QPOs. Although this be- 
havior is complex, the power spectral components 
of the three atoll sources for which we have per- 
formed a multi-Lorentzian timing study show re- 
markable similarities. By using frequency vs. fre- 
quency plots such as the one in Figure 8, and by 
comparing the power spectra of different sources 
directly with each other it is possible to identify 
all these components within a single classification. 
We do this here for 4U 1608-52 (section 3.4). Note 
however, that sometimes additional narrow fea- 
tures become significant in the power spectra (see 
5.2.1). As another example, a recent RXTE obser- 
vation of 4U 1608-52 during an outburst in 2001, 
shows a power spectrum (Figure 15) very simi- 
lar to that of interval B (Fig. 7), except that in- 
stead of the broad peak at ~ 25 Hz, a narrow 
[Q = 5.5) QPO with an rms fractional ampli- 
tude of 4.3 ± 0.6% appears at 30 Hz. Also several, 
marginally significant, features are visible between 
2 and 5 Hz. 

If we compare the low-frequency features of the 
atoll sources with those found in the Z sources, in 
a frequency vs. frequency plot we see that only 
the HBO of the Z sources can be unambiguously 
identified with the L/j of the atoll sources. The fre- 
quency vs. Uu relations of the other low-frequency 
features partly overlap in Figure 8 but are not 
identical, and the relation of the characteristic fre- 
quency of the LFN with Vu is not even the same for 
all Z sources. It appears that for some components 
the frequency vs. Vu relations are not universal. 
There might be small physical differences between 
the atoll sources and Z sources (and within the 
Z sources themselves) that affect these relations. 
So, either these features are different phenomena 
in different types of sources, or something else (e.g. 
neutron star mass, magnetic field strength) affects 
the relations. If this is so, this hidden parameter 
does not change much from source to source, or the 
relations do not depend strongly on this parame- 
ter as the relations are similar, as the relations are 
very similar within each source type. 

5.3. QPO models 

Three different characteristic frequencies can be 
used to test the frequency relations predicted by 
the several QPO models; fmax, fo and fdamped- 
J^max is the frequency at which a Lorentzian con- 



tributes most of its variance per log frequency 
(Belloni et al. 2002) and is the one used in this pa- 
per, vq is the centroid frequency of the Lorentzian 
and can be obtained from our fit parameters z/,nax 
and Q as 1^0 = 2t'inaxQ/\/4Q^"-i-T. According to 
Titarchuk (2002) the vo of the Lorentzian is shifted 
with respect to the eigenfrequency of the oscilla- 
tion due to damping. This shift depends on the 
damping rate which can be estimated from the 
width of the Lorentzian (Titarchuk 2002). We can 
calculate this eigenfrequency, i^dampcd, from :^niax 

and Q as Z^damped = t'maxv'l + 1/(4Q^ + !)• In 
this section we use our fit parameter fmax unless 
stated otherwise. 

The transition layer model (TLM Osherovich 
& Titarchuk 1999; Titarchuk & Osherovich 1999) 
associates the lower kilohcrtz QPO with the Kep- 
lerian frequency of the inner disk edge [vk]- Be- 
tween the neutron star and the Keplerian disk 
a transition layer is present. In the TLM the 
upper kilohertz QPO is produced by radial os- 
cillations of a blob thrown out of the transition 
layer into a magnetosphere. This radial eigen- 
mode or hybrid frequency {vhybrid) relates to vk as 
fhyhrid = (t'x "h 4f2^)^/^, where Q. is the rotational 
frequency of the star's magnetosphere near the 
equatorial plane. This implies that Vhybrid should 
always exceed 2Vl. For 4U 1608-52 Osherovich 
& Titarchuk (1999) found that Vl is always larger 
than 300 Hz and thus Vhybrid > 600 Hz should ap- 
ply. For 4U 1608-52 we find that when there is a 
pair of kilohertz QPOs present (intervals E-H) the 
upper kilohertz QPO frequency ranges from 830 to 
1062 Hz. But the 200-682 Hz single kilohertz QPO 
of intervals A-D, which we identify as the upper 
kilohertz QPO based on frequency-color (Mendez 
et al. 1999; van Straaten et al. 2000; Di Salvo et 
al. 2001), frequency-frequency (Fig. 8) and rms- 
frequency (see Mendez et al. 2001, and Fig. 9) 
correlations, fall well below the 600 Hz predicted 
by the TLM model. Note that neither the use of 
fo, nor I'dampod, shifts these frequencies above the 
600 Hz limit. If these single kilohertz QPOs were 
not the upper, but the lower kilohertz QPOs, this 
would lead to the unlikely scenario that the pres- 
ence of an upper kilohertz QPO would lead to a 
completely different power spectrum for a similar 
lower kilohertz QPO frequency (compare, e.g., C 
with E or D with G in Fig. 7). 

The Uh vs. ve relation for neutron stars and 
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BHCs (sec Psaltis ct al. 1999 for the relation in 
the i/Q representation, and Behoni et al. 2002 for 
the relation in the J^max representation) was re- 
cently extended towards lower frequencies by in- 
cluding 17 white dwarf sources (see Mauche 2002; 
Warner & Woudt 2003, and references therein). 
The points of 4U 0614+09, 4U 1728-34, and 4U 
1608-52 fall on the Vh vs. relation; only the 
point from interval B of 4U 1608-52 deviates. The 
Vfi vs. Vf relation can be fitted with a power 
law with an index close to 1 (Psaltis et al. 1999; 
Mauche 2002). The TLM explains this by assum- 
ing that vg is the Keplerian frequency of the in- 
ner disk (see above), and Vh represents the fre- 
quency of magnetoacoustic oscillations, vma, in 
the disk transition layer. The relation between 
Vh and ve is then a result of a global relation be- 
tween uma of the transition layer and uk at the 
adjustment radius (Titarchuk & Wood 2002). For 
each individual source the TLM predicts an index 
for the Uma vs. uk relation that is steeper than 
one and should be studied separately (Titarchuk, 
Bradshaw & Wood 2001). Note, however, that 
the atoll sources cover a large range of the Vh vs. 
Vf relation, as there are points around i/f = 20 
Hz and around ue = 550 Hz. While 4U 1608- 
52 only contributes points around ve = 20 Hz and 
4U 1728-34 only aroimd u/, = 550 Hz, 4U 0614+09 
contributes both. This suggests that the Psaltis et 
al. (1999) relation is not just a global relation be- 
tween different soiirccs as suggested by Titarchuk 
& Wood (2002) but can also be found within indi- 
vidual sources. Note also that although the TLM 
explains the Psaltis et al. (1999) relation as above 
in Titarchuk & Wood (2002), the HBOs of the 
Z sources, which are part of Psaltis et al. (1999) 
relation, are explained differently, namely as the 
vertical eigenmode of a blob rotating with a Ke- 
plerian frequency thrown into the magnetosphere 
(see above) in other TLM papers (e.g. Titarchuk 
2002). 

In the sonic point beat frequency model (SPBFM; 

Miller, Lamb & Psaltis 1998) the upper kilohertz 
QPO represents the Keplerian frequency at the 
inner disk edge. The lower kilohertz QPO then 
arises from a beat of the Keplerian frequency at 
the inner disk edge with the neutron star spin 
frequency. This led to an early prediction that 
Al/, the frequency difference between the upper 
and the lower kilohertz QPO, is equal to the spin 



frequency and should therefore be constant. Ob- 
servations of several sources, among which 4U 
1728-34 (Mendez & van der Klis 1999) and 4U 
1608-52 (Mendez et al. 1998), showed a signif- 
icant decrease in Av when the upper kilohertz 
QPO frequency increased. Further refinements 
of the SPBFM by Lamb & Miller (2001) could 
explain this decrease in Ai/. Neither the use of 
I'max nor fdamped instead of z/Q makes Av con- 
stant. Note that our measurements of the high - 
frequency kilohertz QPOs {v£ > 500 Hz) are not 
the best to use here as the addition of a lot of data 
leads to an artificial broadening of the Lorentzian 
(see §3.4). Therefore here we used the precise 
measurements of these high-frequency kilohertz 
QPOs made by Mendez & van der Klis (1999) for 
4U 1728-34 and by Mendez et al. (1998) for 4U 
1608-52 to calculate fmax and fdamped- For both 
sources Av increases by less than 1.5 % if ^'max 
is used instead of vq and by less than 3.1 % if 
i^damped is used instead of vq. Note that, just as 
for the TLM, the low frequencies we find for the 
upper kilohertz QPO (see above) are a problem for 
the SPBFM. In the SPBFM the lower bound on 
the upper kilohertz QPO frequency is set by the 
Keplerian frequency at the maximal radius where 
the radiation coming from the neutron star sur- 
face can still remove siifficient angular momentum 
from the gas in the disk so that it falls super- 
sonically to the neutron star surface (Miller et al. 
1998). This radius is about 3Rms, where i?„j,, is 
the radius of the marginally stable orbit. So, to 
reach z/„ < 300 Hz the mass of the neutron star 
has to exceed ~ 2.5 Mq. 

The relativistic precession model (RPM; Stella 
fc Vietri 1998, 1999) assumes that the upper kilo- 
hertz QPO represents the Keplerian frequency of 
the inner disk. The lower kilohertz QPO frequency 
represents the periastron precession frequency of 
the accretion disk, which is assumed to contain 
slightly elliptical orbits. This leads to the predic- 
tion that the frequency difference (Av) between 
the two kilohertz QPOs decreases both at low and 
high kilohertz QPO frequencies (see Figure 1 of 
Stella & Vietri 1999). If the broad Le at low fre- 
quencies (intervals A and B of 4U 1608 52 and 
interval 1 of 4U 0614+09 in van Straaten et al. 
2002) is the lower kilohertz QPO (but see §5.2.1), 
then this predicted decrease in Av at low kilohertz 
QPO frequencies is observed in 4U 0614+09 and 
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4U 1608-52. Note that the point of interval B of 
4U 1608-52 falls above the curves of Figure 1 in 
Stella & Vietri (1999) by 3a. More complicated 
modelling within the framework of the RPM as 
is done for Sco X-1 in Stella & Vietri (1999) is 
beyond the scopes of this paper. 

The RPM also predicts a low-frequency QPO 
at the Lense-Thirring precession frequency of 
ulj, = 4.4 X IQ-^h^M-'^Vsvl. Hz (Stella & Vietri 
1998), where J45 is the moment of inertia in units 
of 10^^ g cm^, M is the mass of the neutron star 
in Mq, and Vs is the neutron star spin frequency. 
We searched our data for a power law with an 
index of 2 as predicted by the RPM by fitting the 
frequency vs. Vu relations for the low frequency 
features (Fig. 8) for the three sources 4U 1728 34, 
4U 0614+09 and 4U 1608 52 with power laws. We 
used Vmayi^ ^0, and i^dampod for the characteristic 
frequencies (see above). We find that the Vh vs. 
Vu relation has an index that is the closest to the 2 
predicted by the RPM. If we use Uraax or i^damped 
for the frequency the relations clearly deviate from 
a power law. The points fall below the fitted power 
law at low and high frequencies and above at inter- 
mediate frequencies, leading to very large x^/dof 
of 4605/21 [v^^) and 7179/21 (i^damped)- If we 
use for the frequency we find that the Vh vs. 
Vu relation can be described by a power law of the 
form V = 5.4(±0.4) x lO-^i/^"^^"^" "^^ (see Figure 
16). The xVdof of this fit is still high (373/20) 
but this is not due to systematic deviations from 
a power law. The averaged relative scatter around 
the fitted power law is only 6 % on the Vh axis 
and only 3 % on the axis. The high x^/dof is 
probably a result of not including any systematic 
errors when determining the errors of the power 
spectral fit parameters (see Ford & van der Klis 
1998). Note, that if we use uq for the frequency we 
exclude those points that have vq = Q (from fits 
with a zero-centered Lorentzian) . The result that 
the Vh vs. Uu relation can be described by a power 
law with an index close to 2 was found previously 
for a more limited frequency range (;/„ > 400 Hz) 
in 4U 1728-34 (index 2.11; Ford & van der Klis 
1998) and 4U 0614+09 (index 2.46; van Straaten 
et al. 2000). It is remarkable that now that we 
have extended the Vh vs. Vu relation towards 
Vu ~ 200 Hz, the relation can still be described by 
a power law and that the power law index is found 
to be even closer to 2. We can get an indication of 



the spin frequency by looking at the burst oscil- 
lation frequency which is 363 Hz for 4U 1728-34 
(Strohmayer et al. 1996) and 620 Hz for 4U 1608- 
52 (Muno et al. 2002b). If these two frequencies 
are the spin frequency of the neutron star in each 
case (the most likely scenario, see Chakrabarty et 
al. 2003), then it is difficult to understand why 
the Vh vs. Vu relations of these sources coincide. 
If we assume that we see the spin frequency in 4U 
1728-34 and twice the spin frequency in 4U 1608- 
52, we find an average spin frequency of about 335 
Hz for these sources. For the vq representation fit, 
this leads to I45/M = 3.7 which is too large for 
proposed equations of state (the acceptable range 
is 0.5 < 745 /M < 2; Stella & Vietri 1998). If 
we assume that we see the spin frequency in 4U 
1608-52 and half the spin frequency in 4U 1728- 
34, the average spin frequency is about 670 Hz for 
these sources. This leads to I45/M = 1.8 which is 
acceptable. 

5.4. The energy spectrum of 4U 1608-52 

Recently Gierlinski & Done (2002b) studied 
the X-ray energy spectra of 4U 1608-52 as a func- 
tion of position in the color-color diagram. They 
find that, similar to other atoll sources, the spec- 
tra of the island state (our intervals A, B and C) 
of 4U 1608-52 are dominated by a hard power law 
spectrum. There is also a soft component present 
which in their view is probably due to the neutron 
star surface rather than the accretion disk. The 
banana state (our intervals E-J) spectrum is soft 
and it is uncertain whether the soft component is 
due to the neutron star surface or the disk. To 
explain the spectral behavior of 4U 1608- 52 Gier- 
linski & Done (2002b) use a scenario where in the 
island state the disk is far from the surface of the 
neutron star. The large radius of the inner disk 
prevents the disk from being observed directly in 
the PCA spectra. In the island/banana state tran- 
sition the inner accretion flow, which was geomet- 
rically thick and optically thin in the island state, 
collapses into a geometrically thin and optically 
thick disk. The inner edge of the disk is now close 
to the neutron star, and it is observed as the soft 
component and increases the soft color. From the 
DISKBB model, Gierhnski & Done (2002b) esti- 
mate an inner disk radius in the banana state of 
^ 30 km. The inner accretion flow becomes more 
optically thick lowering the hard color. We can 
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compare this scenario with the results of our tim- 
ing study. If we assume that Vu represents the Ke- 
plerian frequency at the inner disk edge (e.g. Miller 
et al. 1998) wc find, assuming a neutron star mass 
of 1.4 Mq (for a different neutron mass the results 
change by a factor (M]vs/1.4)^/^, which is 1.2 for 
a 2.4 M0 neutron star), that the inner disk termi- 
nates at 47 — 28 km for the island state (intervals 
A, B and C). For the banana state (intervals E-J) 
we find an inner disk radius of ^ 19 — 15 km and 
decreasing towards higher soft color. Indeed we 
find that the inner disk radius is closer to the neu- 
tron star in the banana than in the island state, 
but the values in both the banana and the island 
state are close to the ~ 30 km estimated by Gier- 
linski & Done (2002b) for the banana state. So, 
the difference in radius between island and banana 
state seems too small to make the disk disappear 
from the PCA spectrum in the island state as pro- 
posed in the scenario of Gierlinski & Done (2002b) 
described above. If we assume that instead of Vu, 
Vf represents the Keplerian frequency at the in- 
ner disk edge (e.g. Titarchuk & Osherovich 1999), 
this would lead to much larger differences in inner 
disk radius between the island (280 — 201 km) and 
the banana state (26 — 20 km), so this interpre- 
tation could be consistent with the Gierlinski & 
Done (2002b) scenario. The identification of as 
the Keplerian frequency at the inner disk edge is 
also suggested by the recently proposed extension 
of the PBK relation (Psaltis ct al. 1999; Belloni 
et al. 2002) towards lower frequencies by includ- 
ing 17 white dwarf sources (see Warner & Woudt 
2003, and references therein). For the white dwarf 
sources a QPO is plotted versus the dwarf nova 
oscillation, DNO, the DNOs are thought to occur 
at the Keplerian frequency of the inner disk edge 
(Warner & Woudt 2002, and references therein). 

The behavior in the soft color vs. intensity di- 
agram of 4U 1608 52 and 4U 1705 44 in their ex- 
treme island states can also be explained within 
the scenario of Gierlinski & Done (2002b). Let us 
assume that represents the Keplerian frequency 
of the inner disk edge. In the extreme island states 
the disk is far out at 200 km and can not be seen 
in the PCA energy spectrum (sec above), the soft 
color should therefore only depend on the spec- 
tral properties of the neutron star. If we look at 
the extreme island states A and B of 4U 1608-52 
and 4U 1705-44 in the soft color versus intensity 



diagram, there seems to be a one to one relation 
between soft color and intensity (see Fig. 6 for 4U 
1608-52). We can fit this relation with a power 
law with an index of about 0.1. We propose that 
in the extreme island state as the accretion rate 
(intensity) increases, the neutron star surface gets 
more heated and the soft color increases. Outside 
the extreme island states the soft color is deter- 
mined by the spectral properties of both the neu- 
tron star surface and the disk and its dependence 
on accretion rate is therefore more complicated. 

6. CONCLUSIONS 

The main conclusions of this paper can be sum- 
marized as follows: 

1 The timing behavior of 4U 1608 52 is al- 
most identical to that of the atoll sources 
4U 1728-34 and 4U 0614-1-09, and can be de- 
scribed in terms of a well defined set of com- 
ponents (L„, Li, LhHz, L/i, Lb, Lb2) which 
all vary in characteristic frequency together 
(except for LhHz for which Vm^x is constant), 
sharing the same frequency correlations. 

2 The timing behavior is not consistent with 
the idea that 4U 1608-52 traces out a three- 
branched Z shape in the color-color diagram 
along which the timing properties vary grad- 
ually as in Z sources. Instead, the power 
spectral properties change smoothly along 
an "e-shaped" track, which is probably com- 
posed of three, partially observed, extended 
extreme island branches. 

3 4U 1608-52, at its lowest frequencies, shows 
a power spectrum almost identical to that of 
the low luminosity bursters IE 1724-3045, 
GS 1826-24, and SLX 1735-269, and the 
milfisecond X ray pulsar SAX J1808.4-3658. 
The low luminosity bursters and the millisec- 
ond X-ray pulsar all seem to be atoll sources 
at low mass accretion rate. A similar con- 
clusion was reached previously based on 4U 
0614-h09 (van Straaten et al. 2002). 

4 The high frequency peak at about 200 Hz 
in interval A of 4U 1608-52 is more likely to 
be the upper kilohertz QPO than the hecto- 
hertz Lorentzian. This means that the corre- 
sponding ~ 200 Hz peaks in 4U 0614-1-09, IE 
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1724 3045, GS 1826 24, SLX 1735 269 and 
SAX J1808.4-3658 also most likely represent 
the upper kilohertz QPO. 

5 The low-frequency part of the power spectra 
of both the atoll and the Z sources behaves 
in a very complex manner. Nevertheless, for 
the atoll sources it is possible to identify all 
these components within a single classifica- 
tion. If we compare the low-frequency fea- 
tures of the atoll sources with those found 
in the Z sources we find that the HBO of 
the Z sources can be identified with of 
the atoll sources. Based upon the frequency 
relations with the upper kilohertz QPO nei- 
ther the LFN nor the sub-harmonic of the 
HBO of the Z sources can straightforwardly 
be identified with the band-limited noise of 
the atoll sources. 

6 The parallel tracks in the frequency versus 
intensity diagram can be identified in the 
intensity versus color diagrams. However, 
the statistical spread in the colors, which 
shows up as the narrow nearly vertical par- 
allel tracks in the color intensity diagrams, 
cause the i^— color correlation to be veiled on 
short timcscalcs (less than hours) so in the 
color intensity diagrams parallel tracks sim- 
ilar to the parallel QPO tracks can only be 
identified on longer timcscalcs. 

7 We have tested the transition layer model, 
the sonic point beat frequency model, and 
the relativistic precession model using our 
results for the three atoll sources 4U 0614+09, 
4U 1728-34, and 4U 1608-52. Neither the 
transition layer model nor the sonic point 
beat frequency model can explain the upper 
kilohertz QPO frequency range we find in 
these sources. The i^h vs. f„ can be de- 
scribed by a power law with an index of 2 
as predicted by the relativistic precession 
model, however it is very likely that the 
burst oscillation frequencies found for 4U 
1728-34 (363 Hz), and 4U 1608-52 (620 Hz) 
represent the spin frequency of the neutron 
star in each case. Then it is difficult to un- 
derstand why the fh vs. f„ relations of these 
sources coincide. 

8 Our timing and color results are consistent 



with the scenario proposed by Gierlinski & 
Done (2002b) to explain the results of their 
X-ray energy spectral study of 4U 1608-52, 
only if and not i/^ represents the Keple- 
rian frequency of the inner disk edge. In this 
scenario the disk is close to the neutron star 
in the banana state, and far from the neu- 
tron star in the island state. The large radius 
of the inner disk in the island state prevents 
the disk from being observed directly in the 
PGA spectra. 
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A. DETAILED DESCRIPTION OF THE POWER SPECTRA 



Wc now provide a clctailecl description of tlie components found in the combined power spectra of 4U 
1608-52 and compare them to previous results on 4U 1728-34, 4U 0614+09 of van Straaten et al (2002) and 
the low luminosity bursters discussed by Belloni et al. (2002). The combined power spectrum of interval A 
(Fig. 7) is very similar to interval 1 of 4U 0614+09 (see Figure 2 in van Straaten et al. 2002). It shows L5 and 
Lft, at the lowest characteristic frequencies observed in atoll sources. Figure 8 confirms the identification of 
these components. As was the case in 4U 0614+09, interval A shows at ~ 15 Hz and a broad Lorentzian at 
~ 200 Hz. The Lorentzian at - 200 Hz found in 4U 1608 52 (interval A) and 4U 0614+09 can be identified 
cLS either L^^ or LhHz based on its frequency. However, this Lorentzian has a similar rms fractional amplitude 
for both 4U 1608-52 and 4U 0614+09. The upper kilohertz QPOs of 4U 1608-52 and 4U 0614+09 always 
have a similar fractional rms, whereas for I^hHz the fractional rms is always about 5 percentage points less 
in 4U 1608-52 than in 4U 0614+09 (see Fig. 9). This hints towards this component being L„. Interval 
B is similar to interval A, it shows L^, L/j, L^ and L„, but all at higher characteristic frequencies. Interval 
B also shows a narrow QPO at 2.458 Hz which will be discussed in §5.2.1. Intervals A and B both show a 
component designated L^, at ~ 15 and ~ 25 Hz respectively. The broad L^ at these low frequencies can be 
tentatively identified with the lower kilohertz QPO based on extrapolations of frequency-frequency relations 
(sec Psaltis et al. 1999; Belloni et al. 2002; van Straaten et al. 2002). However, this Lorentzian might also 
be an unrelated feature (see §5.2.1). Interval C shows only three components, Lf,, and L„ again all 
at higher frequencies than in the previous interval. Interval D shows Lb, LhHz and L„. Compared to 
interval C, is no longer present and h^Hz has appeared. Based on the similarities with 4U 1728 34 and 
4U 0614+09, Lh would be expected to be more narrow than in intervals A-C (expected Q about 1-2) and 
have a characteristic frequency of about 30 Hz. Figure 7 shows a non-significant narrow peak at ~ 21 Hz 
in interval D which can probably be attributed to h^. As noted above, I^hHz, a- broad Lorentzian with a 
constant i^max of ~ 150 Hz, appears in interval D for the first time. In interval A h^Hz would lie on top of 
Lu (see above) and in intervals B and C Lhhz would lie in the flank of L„. Therefore it is hard to detect 
LhHz in intervals A C. Note that there is some excess power at the expected location of 150 Hz in interval 
C. In interval D L„ has moved above this frequency range so that LhHz is more easily detectable. Both the 
frequency and the fractional rms of Lf, in interval D are above the relations in Figures 8 and 9. Including 
the little peak at ~ 21 Hz in the fit hardly affects Li,. These deviations from the relations in Figures 8 
and 9 are most likely due to frequency shifts inside interval D, which covers a large range in the color-color 
diagram (see Fig. 5). Splitting this interval further up would compromise the statistics. Interval E shows 
Lft, Li,2 and a pair of kilohertz QPOs. Here the character of L5 has changed from a band-limited noise 
component (fitted with a zero-centered Lorentzian) to a QPO. LhHz is visible in the power spectrum but is 
not significant (2.4 a single trial) in this interval. Because leaving it out of the fit makes Lb2 fit power up 
to hundreds of Hz leading to an overestimated i^b2 and rms fractional amplitude, we include a LhHz at 150 
Hz (the average value in 4U 1728-34 and 4U 0614+09) with Q = 0.2 (see above) in the fit. In Table 3 and 
in Figure 9 we only show the upper limit for this LhHz- The characteristic frequency of Lf,2 in interval E 
is still high compared to the general relation between iyh2 and ly^ (sec Fig. 8). This is similar to what was 
observed in 4U 1728-34 where in this range there are also two points where i'(,2 is significantly above the 
relation. These deviations all occur close to the Lb "transformation" and when the Lb QPO is still weak 
compared to the band limited noise component Lb2 (see Fig. 9). so one possible explanation is that in these 
transition power spectra a narrow and a broad Lb are merged together. The power spectra of intervals F, 
G and H all show very low-frequency noise (VLFN), and Lf,, Lb2, LhHz and the pair of kilohertz QPOs 
all increasing in frequency. Intervals I and J both show VLFN, LhHz and a broad kilohertz QPO. It is 
unclear whether the kilohertz peak represents the upper, the lower or even a blend of both kilohertz peaks. 
There were no other power spectral features present which could be used to identify this component with the 
help of the correlations of Figure 8. If these peaks, with characteristic frequencies of about 1050 and 1250 
Hz, would be the lower kilohertz QPO they would be among the highest frequency lower kilohertz QPOs 
found in any source up to date (see e.g. van der Klis 2000). So to be slightly conservative and also because 
it is more convenient, as in Figures 8 and 9 we use the characteristic frequency of L„ to plot against, we list 
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these kilohertz QPOs as L„. The Q values for the kilohertz QPOs in intervals I and J are much lower than 
those in 4U 1728-34 and 4U 0614+09. This is probably an artificial broadening due to the movement of the 
kilohertz QPO during the intervals. A finer subdivision of intervals I and J does not lead to observing L„ 
and separately. Note that no kilohertz QPOs were found in the individual observations of these intervals 
(Mendez et al. 1998). Interval J shows an additional broad bump which is also seen in 4U 1728-34 (interval 
18 and 19 of Figure 1 in van Straaten et al. 2002) and 4U 0614+09 (interval 9 of Figure 2 in van Straaten 
et al. 2002); this Lorentzian can be identified as either Lb, L(,2 or a new type of Lorentzian. Here we list it 
as Lfe2. 

The VLFN in 4U 1608-52 shows a similar behavior to that in other atoll sources. The power law index 
is about constant around 1.5, with the exception of interval F where it is 2.5. The rms fractional amplitude 
generally increases with position in the color-color diagram. In intervals I and J the VLFN clearly deviates 
from just a power law. This leads to high x^/dof values of 2.5 and 1.8. Adding Lorentzians in interval I 
and J to fit the deviations leads to significantly better x^/dof of 1.3 and 1.2. The Vma.x of the Lorentzian in 
interval I is then 0.02 Hz and in interval J 0.3 Hz In interval I the Lorentzian does not only fit the deviation 
but almost the entire VLFN. The inclusion of the Lorentzians lowers the rms fractional amplitude of the 
power law to 1.59 % for interval I and 2.27 % for interval J. 
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Table 1 



Interval 


Duration 


At 


Intensity 


Classification 


Number 


(days) 


(days) 


(Crab) 




segment 1 of the lightcurve 


1 


3 


3 


0.16-0.25 


H 


2 


1 


3 


0.05 


D 


3 


1 


3 


0.04 0.05 


I 


4 


1 


3 


0.02-0.03 


D 


5 


1 


4 


0.02 


C 


6 


1 


16 


0.03-0.04 


D 


7 


116^ 


84 


0.03-0.04 


B 


8 


64^ 


402 


0.01-0.02 


C 


segment 2 of the hghtcurve 


9 


43 


5 


0.24-1.20 


J 


10 


6 


1 


0.08-0.21 


FGH 


11 


1 


1 


0.05 


E 


12 


2 


1 


0.02-0.04 


C 


13 


1 


1 


0.03 


E 


14 


5 


2 


0.05-0.10 


I 


15 


1 


7 


0.03-0.04 


FG 


16 


7 


38 


0.04-0.05 


B 


17 


119^ 


523 


0.03-0.05 


A 


segment 3 of the hghtcurve 


18 


1 


4 


0.06-0.07 


I 


19 


1 


2 


0.03 


D 


20 


1 


7 


0.01 


C 


21 


1 


3 


0.03 


J 


22 


1 


9 


0.01 


C 


23 


1 


12 


0.02 


I 


24 


1 


4 


0.01 


C 


25 


1 


4 


0.02-0.03 


G 


26 


1 


4 


0.02 


C 


27 


1 


4 


0.07-0.08 


F 


28 


6 


2 


0.03-0.11 


I 


29 


1 


2 


0.01 


D 


30 


1 




0.01 


C 



^ Large time gaps occured between observations (see Figs. 2 
and 3). 



Note. — For each continuous time interval as defined in §3 we 
list the duration in days (with a minimum value of 1 day), the 
number of days until the next interval (At), the intensity range 
in the 2-16 keV energy band, and a classification based on the 
power spectrum and the position of the source in the color-color 
diagram (see §3.4). For the duration, At and the intensity range 
the 256 s lightcurve obtained in §2 was used. 
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Table 2 

Low-frequency fit parameters of the multi-Lorentzian fit for 4U 1608-52 



Interval 


VLFN 




H 






^62 










Letter 


rms 
(%) 




(Hz) 


Q 


rms 
(%) 


(Hz) 


Q 


rms 
(%) 


i^max 
(Hz) 


Q 


rms 

(%) 


A 






0.247±0.012 


(fixed) 


11.03±0.36 








1,613±0,098 


n.l56±0.091 


12.02±0,91 


B 






1.165±n.045 


(fixed) 


9.97±n.l9 








6.17±0.27 


0.56±0.11 


9.73±0,96 


C 






3.22±0.32 


(fixed) 


9.62±n.67 








19.8±1.5 


0.24±0.14 


11.79±0.86 


D 






22.54±0.94 


(fixed) 


15.09±0.25 












< 6.1^ 


E 






22.77±0.52 




4 46+ "-'J^ 
°-0.64 


30.9±4.9 


(fixed) 


10.4±1.0 






< 2.0^ 


F 


0.840±0.044 


2.49±0.16 


32.11±0.76 


1.15±0.13 


5.99±0.2S 


9.9±2.8 


(fixed) 


3.16±0.42 






< 2.8^ 


G 


0.740±0.078 


1.41±0.16 


37.1±1.1 


1.37±0.22 


4.83±0.31 


10.6+1 = 


(fixed) 


2 44+0-59 
-0.39 








H 


1.180±n.025 


1.537±0.041 


47.4±1.1 


2.0S±0.35 


2.68±0.1S 


i5.4+^:t 


(fixed) 


1.98±0.26 








1 


3.587±0.039 


1.436±0.016 




















J 


2.446±0.019 


1.317±0.ni0 








17.3±1.6 


0.37±0.23 


1.66±0.17 









^95 % confidence upper limits (see §3.4). 



Note. — The fit parameters of the low— frequency part of the power spectra. Listed are the parameters of the very Low-Frequency noise component (VLFN), which is fitted with 
a power law, L^, L52 i and L^ which are fitted with Lorentzians. The horizontal line between intervals D and E marks the "transfer" of Lj, from band— limited noise component 
into a QPO and the simultaneous appearance of L^2 (see §3.4). The horizontal line between intervals I and J marks the broad Lorentzian in interval J for which it is unclear 
whether it is Lj^, Ljj2 °^ ^ new component (see §A). Lj^ and Lj,2 were fitted with a zero— centered Lorentzian (Q = 0) when they showed their band— limited noise like behavior. 
For the power— law a is the power law indices and the fractional rms is integrated from 0.01 to 1 Hz. For the Lorentzians we show the characteristic frequencies (= i-'max)i Q 
values (= fQ/2A) and the integrated fractional rms (from to 00 Hz, over the full PCA energy band). Not listed is the narrow Lorentzian found in interval B, this Lorentzian has 

I'max ^ 2.458 ± 0.037, Q ^ '^■5jl2'l ^"'^ ™^ ^ ^-^^ ^ '^^^ quoted errors use Ax^ ^ i.q. 



Table 3 

High-frequency parameters of the multi-Lorentzian fit for 4U 1608-52 

Interval "Lf^fj^ L^ L^ 



Letter 


I'max 


Q 


rms 




Q 


rms 




Q 


rms 




(Hz) 




(%) 


(h" 




(%) 


"(Szf 




(%) 


A 








14.74±0.82 


0.05±0.10 


15.76±0.87 


216±34 


(fixed) 


14.10±0.58 


B 






< 9.1^ 


24.2±3.9 


(fixed) 


12.48±0.55 


309±13 


0.93±0.14 


12.29±0.52 


C 






< 12. 0'^' 








474 + 21 


0.77+0.14 


18.50+0.73 


D 


189±23 


0.90±0.29 


8.91±0.94 








682.5 + 7.3 


2.96 + 0.27 


16.95 + 0.54 


E 






< 10.3^ 






6.9±1.0 


830.3 + 5.8 


6.11+0.67 


13.64 + 0.54 


F 


134.1±6.1 


1.32±0.27 


5.36±0.34 


602.2±1.3 


19+14 


„ ,+ 1.5 
**-^-0.4 


896.0 + 3.0 


8.30 + 0.61 


8.57 + 0.22 


G 


157±18 


1.17±0.49 
■'■"-0.9 


3.94±0.49 


668.1±2.3 


16.8+4,;" 


9 91 + 0-36 
^■^ -0.20 


973.1 + 6.3 


10.2+1.6 


5.53 + 0.32 


H 


158.9±5.8 


2.13±0.25 


784.35±0.91 


8.04±0.17 


9.549±0.077 


1061.9 + 6.3 


12.5 + 2.8 


3.16 + 0.23 


I 


93.1±5.0 


1.00±0.22 


5.25±0.28 








1046±45l 


2.83±0.75l 


5.98+0.57^ 


J 


123±25 


0.37±0.23 


1.66±0.17 








1252±95l 


0.77±0.22l 


3.61±0.26l 



The identification of this Lorentzian is not certain. 
^95 % confidence upper limits (see §3.4). 

Note. - The fit parameters of the high— frequency part of the power spectra. Listed are the parameters of L^j^^ , , and which are all fitted 

with Lorentzians. The horizontal line separates intervals A and B where L^ is a broad low— frequency Lorentzian (see §3.4). In intervals E— H L^ 
is the lower kilohertz QPO. Note however, that the identification of L^ in intervals A and B as the lower kilohertz QPO is very tentative (see 
§5.2.1). For each Lorentzian we show the characteristic frequencies (= '-'max)) Q values (= l'q/2A) and the integrated fractional rms (from to 
oo Hz, over the full PCA energy band). The quoted errors use Ax^ = 1-0. 
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CD 

c\i 
o ^ 

CV3 
CD 

cvi 

M 

in CV2 
cn 

6 ^ 

in ^ CV2 



O <^ 

8 to 

t; cvi 

o ^ 



CD 

cvi 

2 CV2 
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Fig. 1. — Soft color of the Crab versus time for each PCU. Soft color is defined as the ratio of count rates 
in the energy bands 3.5—6.0/2.0—3.5 keV. Dashed lines indicate the gain changes between gain epochs 1-5. 
The soft color (~ 1.8) in epoch 5 for PCUO, after its propane layer was lost, is offscale. 
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Fig. 2. — Lightcurvc (top) and color-color diagram (bottom) for the decay of the 1998 outburst. In the bottom frame, the 
black points correspond to the 1998 outburst and the grey points to the full data set. Intensity is defined as the count rate in 
the energy band 2.0—16.0 kcV; hard and soft colors arc the count rate ratios in the 9.7—16.0/6.0—9.7 keV and 3.5—6.0/2.0—3.5 
keV bands respectively. Intensity and colors are in units of Crab (sec §2). Each point in the color-color diagram spans 256 s. 
For clarity no errors arc shown, but relative errors are always below 5% (see §2). Each point in the lightcurvc corresponds to a 
one-day average. The numbers indicated in the plot correspond to continuous time intervals for which the power spectra and 
the position in the color-color diagram remain similar. The symbols indicate positions in the color-color diagram for which the 
power spectra are similar; a legend linking the symbols to the classification of §3.4 is provided in the top frame. The enlarged 
crosses at the peak of the outburst and the corresponding enlarged crosses in interval 9 of the color-color diagram indicate data 
for which no power spectra could be computed due to data overflow (see §2). 
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Fig. 3. — As Figure 2, but for the 1996 outburst. 
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Fig. 4. — As Figure 2, but for the third segment of the data (persistent data). In this case, the classification 
marked by the symbols is based on the position in the color-color diagram only (see §3.3). 
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Fig. 5. — Color-color diagram divided into 10 representative intervals marked from A to J. Colors are in 
units of Crab (see §2). 
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Fig. 6. — Soft color (top panel) and hard color (bottom panel) plotted versus intensity. The symbols 
correspond to the 10 representative intervals of Figuryj^S (see §3.4). Intensity and colors are in units of Crab 
(see §2). 




Frequency (Hz) Frequency (Hz) 

Fig. 7. — Power spectra and fit functions in the power spectral density times frequency representation (see 
§2) for 4U 1608-52. The different lines mark the individual Lorentzian components of the fit. The dashed 
lines mark both Lb and L{,2, the dotted lines L/j, the dash-dotted line LhHz, and the dash-dot-dot-dotted 
line Le and L„. In intervals F to J a power-law is included to fit the VLFN; this power-law is indicated 
with a solid line. A narrow Lorentzian found in interval B is indicated with a solid line. The power spectral 
components are also indicated in the plots. Interval letters are indicated. 
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Fiff. 8 Correlations between the eharactcrislic frequencies (— fmax) of the several Lorentzians used to fit the power spectra of 
4U 1608—52 and the z^max of the Lorentzian identified as the upper kilohcrtz QPO, plotted together with the results for 4U 1728-34 
and 4U 0614+09 of van Straatcn ct al. (2002). The black points mark the results for 4U 1608-52, the grey points the results for 4U 
1728—34 and 4U 0614+09. The solid dots mark Lb2. the triangles L;,, the crosses L^, the squares I-ihhz and the diamonds L^. The two 
points with the big circles on the very left are from 4U 0614+09 and interval A of 4U 1608-52, for which the fourth Lorentzian can be 
identified as either the upper kilohertz QPO or the hectohertz Lorentzian based on its frequency (but probably is the upper kilohertz 
QPO; see §A). We use the parameters of this Lorentzian both for the upper kilohertz QPO and for the hectohertz Lorentzian. 
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Fig. 9 . The fractional rms amplitude of all components plotted versus the fmax of the upper kilohertz QPO. The black triangles 

mark our results for 4U 1608—52, the grey stars the results from 4U 0614+09, and the grey squares the results from 4U 1728—34 (van 
Straaten et al. 2002). Arrows indicate 95 % confidence upper limits (see §3.4). The two points with the big circles on the very left are 
from 4U 0614+09 and interval A of 4U 1608—52, for which the fourth Lorentzian can be identified as either the upper kilohertz QPO or 
the hectohertz Lorentzian based on its frequency (but see §A). We use the parameters of this Lorentzian both for the upper kilohertz 
QPO and for the hectohertz Lorentzian. The dashed line in the second panel marks the point below which L^ is a broad low— frequency 
Lorentzian (see §3.4). The dashed line in the fifth panel marks the "transformation" of L^ from band— limited noise component into 
a QPO and the simultaneous appearance of Lf,2 (see §3.4). The dashed line in the sixth panel marks the point above which a broad 
Lorentzian is present for which it is unclear whether it is L^, Lb2 or a new component (see §A). 
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Fig. 10. — The frequency of the lower kilohertz QPO (top panel), hard color (bottom panel) and soft color 
(middle panel) plotted versus 2.0-16.0 keV intensity. Intensity and colors are in units of Crab (sec §2). 
The kilohertz QPO frequencies are from Mcndcz ct al. (2001). Similarly to what was done in Mendez et 
al. (2001) we only include data for which both kilohertz QPOs are detected simultaneously and therefore 
the lower kilohertz QPO is identified unambiguously. The symbols and the alternating black/grey represent 
the parallel tracks in the intensity versus frequency diagram that are continuous in time and only have gaps 
due to earth occultations of the source (~ 2500 s). Typical soft and hard color error bars are shown in the 
middle and bottom panel. The errors in the intensity are smaller then the size of the symbols. 



37 




0.95 1 1.05 

Soft color (Crab) 

Fig. 11. — The frequency of the lower kilohertz QPO plotted versus soft color. The kilohertz QPO frequencies 
are from Mendez et al. (2001). Lines connect data from each individual satellite orbit. For clarity four 
representative individual orbits are highlighted in black. Soft color is in units of Crab (see §2). A typical 
soft color error is shown. 
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Fig. 13. — Characteristic frequency of \jh versus that of the narrow low-frequency QPO. The dashed hne 
indicates a power law fit to the low-luminosity bursters IE 1724-3045 and GS 1826-24, and the BHC GX 
339-4. The dash-dotted line is a power law with a normalization half of that of the dashed line. 
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Fig. 14. — Characteristic frequencies (see §5 for definitions) of the different power spectral components 
versus i^max of the upper kilohertz QPO for atoll and Z sources. The grey symbols are the results of the 
multi-Lorentzian fit to the atoll sources 4U 1608 52, 4U 0614+09, and 4U 1728-34 of Figure 8. The black 
symbols represent the Z sources; GX 5-1 (Jonker et al. 2002), GX 340+0 (Jonker et al. 2000), GX 17+2 
(Homan et al. 2002), and Cyg X-2 (Kuznetsov 2002). For the Z sources we include, as indicated in the plot, 
the lower kilohertz QPO, the Low-Frequency Noise (LFN), the Horizontal Branch Oscillations (HBO), and 
the harmonic and sub-harmonic of the HBO. 
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Fig. 15. — Power spectrum of 2001 November 23 in the power spectral density times frequency representation 

(sec §2) for 4U 1608 52. This power spectrum is very similar to that of interval B (see Figure 7), except 
that instead of the broad peak at ~ 25 Hz, a narrow {Q = 5.5) QPO appears at 30 Hz. 
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